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Experimental preparation of eight-partite linear and two-diamond shape cluster states 

for photonic qumodes 

Xiaolong Su, Yaping Zhao, Shuhong Hao, Xiaojun Jia, Changde Xie, and Kunchi Pen£0 

State Key Laboratory of Quantum Optics and Quantum Optics Devices, 
Institute of Opto-Electronics, Shanxi University, Taiyuan, 030006, People's Republic of China 

The preparation of multipartite entangled states is the prerequisite for exploring quantum in- 
formation networks and quantum computation. In this letter, we present the first experimental 
demonstration of eight-partite spatially separated CV entangled states. The initial resource quan- 
tum states are eight squeezed states of light, through the linearly optical transformation of which two 
types of the eight-partite cluster entangled states are prepared, respectively. The generated eight 
entangled photonic qumodes are spatially separated, which provide valuable quantum resources to 
implement more complicated quantum information task. 

PACS numbers: 03.67.Bg, 03.67.Lx, 03.65.Ud, 42.50.Dv 



Developing quantum information (QI) science have exhibited unusual potentiality 0, Optical QI based on 
^ \ exploiting discrete-variable (DV) of single-photon states (photonic qubits) and continuous-variable (CV) of optical 
modes (photonic qumodes) plays important role in QI development. The one-way quantum computation(QC) based 
on multipartite cluster entanglement is initially proposed by Raussendorf and Briegel in the DV model [3], then it is 
extended to the CV regime by Menicucci et al [4J. For one-way QC model the qubits (qumodes) are initialized in 
a multipartite cluster entangled state firstly, then a variety of quantum logical operations can be achieved only via 
O |. the single-qubit (qumode) projective measurement and the classical feedforward of the measured outcomes, in which 
! 1 1 ' the order and choices of measurements are determined by the required algorithm Q. The basic logical operations of 
£h ! one-way DVQC has been experimentally demonstrated by several groups @t3|- 

Parallelly, the theoretical and experimental explorations on one-way CVQC were also proceeding continually [H- 
Tol. fl^4l5j . In contrast of the probabilistic generation of photonic qubits in most cases, CV cluster states are produced 
in an unconditional fashion and thus the one-way QC with CV cluster entangled photonic qumodes can be imple- 
mented deterministically p^4l9j . Following the theoretical proposals on one-way CVQC the principally experimental 
demonstrations of various one-way QC logical operations over CVs were achieved by utilizing bipartite and four-partite 
cluster entangled photonic qumodes, respectively fl2l - fl5j ]. To develop more complicated QC larger cluster states with 
more numbers of entangled qubits (qumodes) are desired. However, the numbers of spatially separable entangled 
qumodes generated by experiments still stay below four-partites, so far [l6l~18j. In the paper, we present the first 
experimental achievement on producing CV eight-partite entangled states for photonic qumodes. Using eight squeezed 
• ■ states of light to be the initial resource quantum states and passing through the linearly optical transformation on a 
specially designed beam-splitter network, the eight-partite linear and two-diamond shape cluster states for photonic 
. qumodes are prepared, respectively. The entanglement feature among the obtained eight space-separated photonic 
' qumodes is confirmed by the fully inseparability criteria of CV multipartite entangled states proposed by van Loock 
£J \ and Furusawa [J]. 

. £^ ■ The cluster state is a type of multipartite quantum entangled graph states corresponding to some mathematic 
graphs 0, 0, E3] • The CV cluster quadrature correlations (so-called nullifiers) can be expressed by 0, [H, Ell 



{p a - V x b ) -> 0, V a£G (1) 



where x a = (a + oJ)/2 and p a — (a — a* ! )/2i stand for quadrature-amplitude and quadrature-phase operators of an 
optical mode a, respectively. The subscript a (b) expresses the designated mode a (b). The modes of a £ G denote 
the vertices of the graph G, while the modes of b £ N a are the nearest neighbors of mode a. For an ideal cluster state 
the left-hand side of equation (1) trends to zero, which stands for a simultaneous zero eigenstate of the quadrature 
combination [l(|. The CV cluster quantum entanglements generated by experiments are deterministic, but also are 
imperfect, the entanglement features of which have to be verified and quantified by the sufficient conditions for the 
fully inseparability of multipartite CV entanglement [l6l - [l9j . There are different correlation combinations [left-hand 
side of equation (1)] in a variety of CV cluster multipartite entangled states, which reflect the complexity and rich 
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FIG. 1: The graph representation of eight-partite cluster states, a: linear cluster state, b: two-diamond shape cluster state. 
Each cluster node corresponds to an optical mode. The connected lines between neighboring nodes stand for the interaction 
among these nodes. 



usability of these quantum systems. The expressions of the milliners for different graph states depend on their graph 
configurations. 

Figure 1 (a) and (b) show the graph representations of CV eight-partite linear (a) and two-diamond (b) shape 
CV cluster states, respectively, each node of which corresponds to an optical mode and the connection lines between 
neighboring nodes stand for the interaction between the connected two nodes. From equation (1) and Fig. 1, we can 
write the milliners of the linear and the two-diamond shape CV cluster states, respectively, which are pLi — x^ — Sli, 

PL2 - XLl ~ X L3 = S L2 , PL3 ~ %L2 - X L 4 = S L3 , p LA - X L3 - X L5 = 8 L i, PL5 ~ %Li - X L 6 = $L5, PL6 - XL5 - X L7 = S L 6, 

PL7 ~ xls - xls = $L7, PL8 — %L7 — $L8 for the linear states; and p m - x D3 - x Di = S m , PD2 - x D3 - x Di = S D2 , 

PD3 - %Dl - X D2 = S D3 , p04 - X D1 - X D 2 ~ X D5 = 5 D A, PD5 ~ X Di - X D7 - X D $ = 5 D5 , p m - X D7 - X D $ = 5 D6 , 

Pm — XD5 — &D6 = $D7, Pds — xd5 — xd6 = ^r>8 f° r the two-diamond states, where the subscripts and 
(i = 1,2, ...8) denote the individual nodes of the linear and the two-diamond shape cluster states, respectively, 5^ 
and express the excess noises resulting from the imperfect quantum correlations. When the variance of 8^ (<5d;) 
is smaller than the corresponding quantum noise limit (QNL) determined by vacuum noises, the correlations among 
the combined optical modes is within the quantum region, otherwise the quantum correlations do not exist. 

The schemes of generating CV multipartite entangled states commonly used in experiments are to achieve a linearly 
optical transformation of input squeezed states on a specific beam-splitter network 1]. Assuming ai and Uki stand for 
the input squeezed states and the unitary matrix of a given beam-splitter network, respectively, the output optical 
modes after the transformation are given by bk = 53/ Uki&i, where the subscripts 1 and k express the designated input 
and output modes, respectively. In our experiment, four quadrature-amplitude i-squeezed states, a m — e~ r im + 
ie +r p"m (m = 1,3,5,7), and four quadrature-phase p-squeezed states, a n — e +r x\? + ie~ r pn ^ (n = 2,4,6,8), are 

applied, where x^ andp^ '' denote the quadrature-amplitude and the quadrature-phase operators of the corresponding 
vacuum field, respectively, r is the squeezing parameter to quantify the squeezing level, r = and r = +oo correspond 
to the two cases of no squeezing and the ideally perfect squeezing, respectively. The unitary matrix for generating 
the CV eight-partite linear cluster state by combining eight squeezed states on optical beam splitters equals to (see 
Supplementary Material) 




The unitary matrix in equation (2) expresses an optical transformation on a beam- 
splitter network consisting of seven beam splitters and can be decomposed into Ul = 
i^C-l^^/st-l^B^l/^ where 
Fk denotes the Fourier transformation of mode k, which corresponds to a 90° rotation in the phase space; B^(Tj) 
stands for the linearly optical transformation on the jth beam-splitter with the transmission of Tj (J = 1, 2, 3 7), 
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where {B^) kk = y/l - T, (B^) kl = y/T, (Bj^) lk = ±VT, and {B±) u = =Fa/1 - T, are elements of beam-splitter 
matrix. /*(— 1) = e™ corresponds to a 180° rotation of mode k in the phase space. 




FIG. 2: Schematic of experimental setup for CV eight-partite cluster state generation. T: transmission efficient of beam splitter, 
Boxes including i are Fourier transforms (90° rotations in phase space), —i is a —90° rotation, and —1 is a 180° rotation, BHD: 
balanced homodyne detector. 

Figure 2 shows the schematic of the experimental set-up for preparing the eight-partite CV linear cluster state. The 
four ^-squeezed and four p-squeezed states are produced by four NOPAs pumped by a common laser source, which 
is a CW intracavity frequency-doubled and frequency-stabilized Nd:YAP/LBO(Nd-doped YAIO3 perorskite/lithium 
triborate) with both outputs of the fundamental and the second- harmonic waves 23] . The output fundamental wave 
at 1080 nm wavelength is used for the injected signals of NOPAs and the local oscillators of the balanced homodyne 
detectors (BHDs), which are applied to measure the quantum fluctuations of the quadrature-amplitude and the 
quadrature-phase of the output optical modes [lf|. The second-harmonic wave at 540 nm wavelength serves as the 
pump held of the four NOPAs, in which through an intracavity frequency-down-conversion process a pair of signal and 
idler modes with the identical frequency at 1080 nm and the orthogonal polarizations are generated pi. l2~ij . Since the 
amplitude and the phase quadratures of the signal and the idler modes are entangled each other, the two coupled modes 
of them at ±45° polarization directions both are the squeezed states @, EH- In our experiment, the four NOPAs are 
operated at the parametric deamplification situation, i.e. the phase difference between the pump field and the injected 
signal is (2n+ 1)tt [n is an integer). Under this condition, the coupled modes at +45° and -45°polarization directions 
are the quadrature-amplitude and the quadrature-phase squeezed states, respectively @,[l6j]. When the transmissions 
of the seven beam splitters are chosen as Tl = 25/34, T2 = T3 = 2/5, T4 = T5 = 1/3, T6 = T7 = 1/2, the 
eight output optical modes bj (j = 1, 2, ...8) are in a eight-partite CV linear cluster state. The quadrature-amplitude 

and quadrature-phase of each bj are measured by eight BHDs, respectively. The nullifiers of the eight output modes 
depend on the squeezing parameters of the resource squeezed states. For our experimental system all four NOPAs have 
the identical configuration (the construction of NOPA is described in the Supplementary Material) and are operated 
under the same conditions. Each of NOPAs is also adjusted to produce two balanced squeezed states. So, the eight 
initial squeezed states own the same squeezed parameter r. In this case we can easily calculate the excess noises of 
the nullifiers for the eight-partite linear CV cluster state consisting of the eight output modes bj (j = 1, 8), which 

are 5 L1 = yf2e^xf\ 5 L2 = y/3e^p { °\ S L3 = -fee^S™ - yfie^ , S Li = ^e-p< 0) + yfle-^ + yj^e^, 

S L5 = Y^e-^i 0) -^e-4 0) -^e-4 0) , S Le = ^"^-^e^, S L7 = -y/te-'x™ and S LS = -V2e~^\ 
respectively. 

The unitary matrix of the two-diamond cluster state Ud equals to UfUl, with Uf = diag{— 1, —i, i, 1, 1, i, —i, — 1} 
(see Supplementary Material), thus the two-diamond shape cluster state can be prepared from the linear cluster state 
via local Fourier transforms and phase rotations. The excess noise terms of the nullifiers of the two-diamond shape 
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cluster state are expressed by 6m = — -^e~ r S:^ + \Zl e_ '^3°'' &D2 = ~~} 2 e ~ r Xi^ ~~ \J\ e ~ r x 3°^ &D3 = ~V^ e ~ r P2 > \ 

<5d7 = \f\ e ~ r jffi ~ ~^2 e ~ T P& ^ anc ^ ^-D8 = \[\ e ~ r ^6°' + "72 e ~ r P8°' ) ' respectively. According to the inseparability 
criteria for CV multipartite entangled states proposed by van Loock and Furusawa [|[ , we deduced the inseparability 
criterion inequalities for CV eight-partite linear and two-diamond shape cluster states, which are given by equations 
(3a)- (3g) and equations (4a)- (4i), respectively (see Supplementary Material). 
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where left-hand sides and right-hand sides of these inequalities are the combination of variances of nullifiers and the 
boundary, respectively. In the Supplementary Material, we numerically calculated the dependencies of the combina- 
tions of the correlation variances in the left-hand sides of equations (3) and equations (4) on the squeezing factor r 
for g = 1 and g = g opt (g opt is the optimal gain factor), respectively. It can be seen, when the optimal gain factors 
are used the correlation variance combinations are always smaller than the boundary for all values of r > 0, i.e. the 
eight-partite CV cluster entanglement always can be realized by the presented system whatever how low the squeezing 
of the initial squeezed states is. However, when g = 1, a lower limitation of r is required to result in the correlation 
variances to be smaller than the boundary. 

The experimentally measured initial squeezing degrees of the output fields from four NOPAs are 4.30 ± 0.07 dB 
below the QNL which corresponds to the squeezing parameter r = 0.50 ± 0.02 (see Supplementary Material). During 
the measurements the pump power of NOPAs at 540 nm wavelength is ~ 180 mW, which is below the oscillation 
threshold of 240 mW, and the intensity of the injected signal at 1080 nm is 10 mW. The phase difference on each 
beam-splitters are locked according to the requirements. The light intensity of the local oscillator in all BHDs is set 
to around 5 mW. The measured QNL is about 20 dB above the electronics noise level, which guarantees that the 
results of the homodyne detections are almost not affected by electronic noises. 

The correlation variances measured experimentally are shown in figure 3 for the linear cluster and figure 4 for 
the two-diamond cluster. They are V{p Ll - x L2 ) = -2.67 ± 0.06 dB, V(p L2 - x L1 - x L3 ) = -2.65 ± 0.13 dB, 
V{p L3 - x L2 - x L4 ) = -2.52 ± 0.20 dB, V{p L4 - x L3 - x L5 ) = -2.69 ± 0.09 dB, V(p L5 - x L4 - x L6 ) = -2.68 ± 0.08 
dB, V[PL6 ~ x L5 - x L7 ) = -2.56 ± 0.10 dB, V(p L7 - x Le - x LS ) = -2.22 ± 0.09 dB, V( PL8 - x L7 ) = -2.21 ± 0.09 
dB and V(pm ~ x D3 - x D4 ) = -2.61 ± 0.10 dB, V(p D2 - x D3 - x D4 ) = -2.57 ± 0.09 dB, V{p D3 - x m ~ x D2 ) = 
-2.39 ± 0.06 dB, V(p D4 - x m - x D2 - x D5 ) = -2.58 ± 0.09 dB, V(p D5 - x D4 - x D7 - x D8 ) = -2.61 ± 0.09 dB, 
V(p D6 - x D7 -x D8 ) = -2.52 ± 0.10 dB, V(p D7 - x D5 - x m ) = -2.59 ± 0.09 dB, V(p D8 - x m - x m ) = -2.58 ± 0.10 
dB, V(p D4 - g D 6Xm - 9D6X D2 - x D5 ) = -1.57 ± 0.09 dB, V{p D5 - x D4 - gDt>x D7 - gDt>x D8 ) = -1.53 ± 0.09 dB. 
From these measured results we can calculated the combinations of the correlation variances in the left-hand sides 
of the inequalities (3a)-(3g) and (4a)-(4i), which are 0.68 ± 0.02, 0.83 ± 0.02, 0.82 ± 0.02, 0.81 ± 0.02, 0.82 ± 0.02, 
0.87 ± 0.02, 0.75 ± 0.02, for the linear cluster and 0.84 ± 0.02, 0.85 ± 0.02, 0.96 ± 0.02, 0.97 ± 0.02, 0.95 ± 0.02, 
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FIG. 3: The measured noise powers of eight-partite linear cluster state. The upper and lower lines in all graphs are shot 
noise level and correlation variances of milliners, respectively, (a)-(h) are noise powers of V(j>li — XL2), V(pL2 — xli — im), 

V{pL3 ~ X L2 - X Li ), V(pLA ~ £ L3 ~ X L5 ), V(pL5 - %L4, - £lb), V(pL6 ~ %L5 - Xlt), V(pL7 ~ %L6 ~ X L s), and V(pLS ~ XL?), 

respectively. The measurement frequency is 2 MHz, resolution bandwidth is 30 kHz, and video bandwidth is 100 Hz. 




FIG. 4: The measured noise powers of eight-partite two-diamond shape cluster state. The upper and lower lines in all graphs 
are shot noise level and correlation variances of milliners, respectively, (a)-(j) are noise powers of V(pm — %D3 — &da), 

V(pD2 — XD3 — XDi), V(pD3 ~ XD1 — Xdz), V(pD4 — Xni ~ X D 2 — Xdb), V(pD5 ~ X D 4 ~ %D7 - XDs), V(pD6 ~ X D 7 ~ &Ds), 

V(pD7 — ins - XD&), V(pD8 — xd5 — x D6 ), V(pD4 ~ ffwifli — ^06^02 — x D5 ), V(pD5 - x D 4 - ffD6*c7 - gDeXDs), respectively. 
The measurement frequency is 2 MHz, resolution bandwidth is 30 kHz, and video bandwidth is 100 Hz. 



0.96 ± 0.02, 0.96 ± 0.02, 0.83 ± 0.02, 0.83 ± 0.02 for the two-diamond cluster, respectively. All these values are smaller 
than the boundary. It means that the prepared two types of CV cluster states satisfy the inseparability criteria for 
verifying multipartite CV entanglement, so the spatially separated eight-partite entangled states of photonic qumodes 
are experimentally obtained. In the experiment we detected the correlation variances under gLi-LS = 1; 9di—D5 = 1 
and goe — 5^6 = 0.60 ± 0.02. For our system, the total transmission efficiency of squeezed beams are about 87% and 
the detection efficiency is about 90%, which lead to the efficient squeezing parameter is r e — 0.30 which is smaller 
than the initially measured squeezing parameter. When the gain factors except goe are taken as 1 and only g°^ is 
utilized, all inequalities in equation (3) and equation (4) are satisfied. If r e > 0.35, the unite gain factor of goe = 1 
can be chosen also (see Supplementary Material). 

In the conclusion, we have experimentally prepared two types of spatially separated eight-partite CV cluster entan- 
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gled states for photonic qumodes by using eight quadrature squeezed states of light and a specifically designed optical 
beam-splitter network. The multipartite entangled states are the essential resources to construct a variety of CVQI 
networks. So far, the single-mode squeezed states over 12.7 dB [13] and the two-mode squeezed states over 8.1 dB 
(28j have been experimentally generated, respectively, based on which and using the presented scheme the CV cluster 
states with more space-separable qumodes and higher entanglement can be obtained. The complexity and versatility 
of CV multipartite entanglement for photonic qumodes not only offer richly potential applications in QC and QI, but 
also provide the basic and handleable entangled quantum states which can be an important tool for further studying 
the amazing and attractive quantum entanglement phenomena. 

This research was supported by the National Basic Research Program of China (Grant No. 2010CB923103) and 
NSFC (Grant Nos. 11174188, 61121064). 
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UNITARY MATRIX FOR GENERATING CV EIGHT-PARTITE CLUSTER STATES 



According to the proposal of Peter van Loock et al [l| , CV cluster states of photonic qumodes can be created via 
a general linear-optics transformation of p-squeezed input modes. If di = e +r x l f" > + ie~ r pf^ and a unitary matrix U 
stand for the annihilation operator of the input modes and the linear-optical transformation respectively, the output 
modes after the transformation are expressed by bk = JTJ; Uki&i, which are the CV cluster states 0, 0. The CV 
cluster states satisfy Ilm[UBi n ] — ARe[UBi n ] — > in the limit of infinite squeezing, where / is the identity matrix, 
B in = ( &i, a,2, ■ ■ ■ , a n ) T is the matrix of input states, A is the adjacency matrix [3j. So we have IlrnU = ARelJ, 
and the unitary matrix is obtained 



U = (I + iA)ReU. 
Based on the unitarity of matrix U, UU^ — I, we have 



(5) 



ReU(ReU) T = (I + A 2 )~ 

In this case, we can obtain ReL 7 and U from the adjacency matrix A. 
For n-partite cluster state, assuming 



(6) 



ReU 



(7) 



where af = ( an, a a , 



ati n ) is a real vectors. According to Eq. (2), we have afdj = (I+A 2 )- 1 = 1, 



where the numbers of these equations are n(n + 1)/2 according to the symmetry of matrix. Since there are n 2 unknown 
numbers in all these equations, we need n(n — l)/2 conditions to solve the equations. For simplicity and without 
lossing generality, some unknown numbers in the equations are chosen to be when we solve the equations. 
For CV eight-partite linear cluster state, the adjacency matrix can be written as 



A = 



(Q 1 \ 

1 1 
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1 1 

1 1 

1 1 
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\ 1 / 



(8) 



we have 
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V 3T 

The matrix elements in Eq. (3) are obtained by the following way. Considering the symmetry, we start from 
the middle row, aj = ( 041, 0:42, 0:43, a 44 , a 45 , a 46 , a 47 , a 48 ). We apply seven initial conditions, 041 = (242 = 
«43 = 0:44 = a 46 = a 47 = Q!48 = 0, then it becomes = ( 0145 ). According to equations a^oti = 

(I + A 2 ) 44 = || , first unknown number a 45 = — ^/§f is obtained. Then we apply six conditions, a 5 i — a 52 — "53 = 

a 56 = a 57 = a 5S = on 5g , we have — ( 0:54 055 ). Using equations afa^ = (I + A 2 )§ 5 = || and 

a^a^ = (I + A 2 )^ = 0, two unknown numbers a 55 = 0, 0154 = — ^/|f are gotten, and so on. 

After all the unknown numbers in Eq. (3) are obtained, the unitary matrix is given from Eq. (1). The unitary 
matrix for eight-partite linear cluster state with eight p-squeezed states to be the input states are expressed 




In the experiment, we prepared four x-squeezed states, a m = e r x m ' + ie +r p m ' (m = 1, 3, 5, 7), and four p-squeezed 
states, a n = e +r x^ + ie~ r pn \n = 2, 4, 6, 8) with four NOPAs, respectively. The transformation between x-squeezed 
state and p-squeezed state can be achieved via a Fourier transformation. Applying Fourier transformation on modes 
o-ii 0-37 &5 and o-7i which corresponds to multiplying i to the values of columns 1, 3, 5 and 7 in the unitary matrix 
we obtain the unitary matrix of CV eight-partite linear cluster state for our experimental system, which is 




Similarly, the unitary matrix for generating CV eight-partite two-diamond shape cluster state can be obtained with 
the same way. The calculated results show that Ud equals to U fUl, where Up = diag{—l, —i, i, 1, 1, i, — i, —1}. The 
unitary matrix Ud is given by 
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II. INSEPARABILITY CRITERIA 



According to the inseparability criteria for CV multipartite entangled states proposed by van Loock and Furusawa 
[H, we deduced the concrete inseparability conditions for CV eight-partite linear cluster state, which are Eqs. (3) in the 
main text. For any separable quantum state, its total density operator can be written as p = JV ViPi,k,...,m <8> Pi,i,...,n 
with a distinct pair of "separable modes" (m, n) and the other modes k ^ I [see equation (25) in Ref. 2.], where 
rji represent the mixture of these separable states. For any combinations u — h\x\ + h%x% + • • • + h^x^ and 
v = gipi + g2Vi + ■ • • + 9nPn, the inseparability criteria are expressed by 0| 



V(u) + V{v) < -( h m g m + ^2 h k9k + h n9n + ^2 higi ) 



(13) 



Based on this equation, we deduced the inseparability criteria for CV eight-partite linear and two-diamond shape 
cluster states, which are Eqs. (3) and (4) in the main text, respectively. 

The violation of each inequality in Eqs. (3) can be used to verify variously partially separable states, that are 



violation of (3a) =>■ p 

violation of (36) =£• p 

violation of (3c) p 

violation of (3d) ==>■ p 

violation of (3e) =>■ p 

violation of (3/) =>■ p 

violation of (3g) =>■ p 



* — 

E, 



Pi,k,. 



(14) 



Pi,k,...,2 69 pi,l 
Pi,k,...,3 <S) pi,l 
Pi,k,...A ® 
Pi,k,...,5 <S) pi,l 
Pi.k 6 ® Pi,Z 



ViPi,k,...,7 59 Pi,/ 



If all inequalities in Eqs. (3) are satisfied, we can confirm that all eight states are inseparable and thus they construct 
a CV eight-partite linear cluster entangled state. 

We also deduced the inseparability criteria for CV eight-partite two-diamond shape linear cluster state, which are 
Eqs. (4) in the main text. Similarly, if the inequalities are violated the partially separable states are expressed by 
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FIG. 5: The dependence of eight-partite linear cluster state to squeezing parameter, a-d are corresponding to inequalities 
(3a)-(3d), respectively. Lines 1 and 3 are left-hand sides of inequalities with unit gain and optimal gain, respectively, line 2 are 
right-hand sides of inequalities. 



violation of (4a) =>■ p 

violation of (46) =£• p 

violation of (4c) p 

violation of (4d) =>■ p 

violation of (4e) =>■ p 

violation of (4/) =>■ p 

violation of (Ag) p 

violation of (Ah) =>■ p 

violation of (Ai) =>■ p 



EViPi,k,...,l <8> pi.l 
i 

22.11iPi,k,-,2 ®Pi,l 

V" ViPi,k,...,l <S> pi,l 
L — 'i 

^ ViPi,k,...,2 <8> pi,l 
* — 'l 

^ rjipi,k,...A ® 

z — *i 

EViPi,k,...,5 ® Pi,; 

Eriipi,k,-,5 ® 
z 

E^ipz.fc 6 ® PiJ 

/^.??iPi,fe,...,6 ® 



(15) 



When all inequalities in Eqs. (4) are satisfied CV eight-partite two-diamond shape cluster entanglement is demon- 
strated. 

Calculating the minimal values of the left-hand sides of the inequalities versus the gain factors (<?li-L8 for the 
linear shape and Qdi-dq for the two-diamond shape) we can obtain the optimal gain factors <?£i-£8 and g°j^l_ D& for 
achieving the detections of the minimal correlation variances. For the linear (<?£i-L8 ) and the two diamond (.o^fi^g) 

_ 13(e 4r -l) _ _ S,(e iT -l) 



cluster states the optimal gain factors are g^i = g^ = ^i^r i 9L2 



9L4 - 

gD6 



9L5 = i 9+ i 5e 4. and g m - w+15eir , g D 2 



_ 15(e 4r 

4(e 4T, - i 
13+4C 



_ 21(e 4r -l) 
— 13+21e 4 '' ' 5-D3 



9(e 4r ~l) 



3L7 - 21 + 13e^ ' 5L3 - 

9(e 8r -l) 
3D4 - 7+18e 4. +9e ^ , gD5 



9L6 — g +8e 4r 
_ 3(3e 8r -2e 4, '-l) 



g+gJB^, i/Z)4 — 7+18e 4 ^+9e 8 '' ' — 7+18e 4 ''+9e s '' ' 

■zr, respectively. The dependence of inseparability criteria of CV eight-partite linear cluster state on the 
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FIG. 6: The dependence of eight-partite two-diamond shape cluster state to squeezing parameter, a, b and c are corresponding 
to inequalities (4a), (4c), and (4e), respectively. Lines 1 and 3 are left-hand sides of inequalities with unit gain and optimal 
gain, respectively, line 2 are right-hand sides of inequalities. 

squeezing parameter is shown in figure 1. In Fig. 1, (a), (b), (c) and (d) correspond to inequalities (3a), (3b), (3c) 
and (3d) , respectively. Lines 1 and 3 are the variance combinations in the left-hand sides of inequalities with unit gain 
and the optimal gain, respectively. Line 2 stands for the boundary corresponding to right-hand sides of inequalities. 
Since the other inequalities have same variances with (3a)- (3c), we only show the four results here. For inequalities 
(3a)-(3d), if the unit gain factor is used only when r > 0.11, r > 0.20, r > 0.24, and r > 0.27 these inequalities 
are satisfied. However, if the optimal gain factors are chosen, these criterion inequalities are always satisfied for any 
r > 0, i.e. the presented system can prepare the CV eight-partite cluster state without the lower limitation of the 
squeezing degree for the input squeezed state. 

The dependence of inseparability criteria of CV eight-partite two-diamond shape cluster state to the squeezing 
parameter is shown in figure 2. Since the variances of inequalities (4a), (4b), (4h) and (4i) are the same, and that of 
(4c), (4d), (4f) and (4g) are also the same, we only plot the correlation variances of the inequalities (4a), (4c) and 
(4e) in Fig. 2a, 2b and 2c, respectively. From Fig. 2a, 2b and 2c, we can see that when the squeezing parameters 
r > 0.20 for Fig. 2a, r > 0.28 for Fig. 2b, and r > 0.35 for Fig. 2c, the correlation variances (curve 1) are smaller 
than the boundary for the case of g = 1. It means that if taking g = 1, a lower limitation for the squeezing parameter 
is required to satisfy each of the inseparability criteria. However, if taking the optimal gain factor g opt (curve 3) all 
variances are below the boundary for any values of r > 0. 

III. EXPERIMENTAL DETAILS 

The four NOPAs were constructed in identical configuration, each of which consists of an a-cut type-II KTP crystal 
and a concave mirror The front face of the KTP was coated to be used for the input coupler and the concave 
mirror serves as the output coupler of the squeezed states, which was mounted on a piezo-electric transducer for 
locking actively the cavity length of NOPA on resonance with the injected signal at 1080 nm. The transmissions of 
the input coupler for the pump laser at 540 nm and the injected signal at 1080 nm are 99.8% and 0.04%, respectively. 
The transmissions of the output coupler at 540 nm and 1080 nm are 0.5% and 5.2%, respectively. The finesses of the 
NOPA for 540 nm and 1080 nm are 3 and 117, respectively. Through the intracavity parametric down conversion 
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FIG. 7: The measured noise power of initial squeezed state, a and b are amplitude and phase quadrature, respectively. Red 
and black lines are the noise power for the QNL and the squeezed state, respectively. 



process of type-II phase match, the two quadrature squeezed states of light at 1080 nm were produced [6j, |7j. The 
squeezing parameter r depends on the strength and the time of parametric interaction in the NOPA. In the calculations 
we have assumed that squeezing parameter for the eight squeezed states is identical for simplicity and without lossing 
the generality. The requirement is easy to be reached in the experiments if the four NOPAs were constructed in 
identical configuration and the intracavity losses of the signal and the idler modes in each NOPA were balanced. The 
experimentally measured squeezing degrees of the output fields from four NOPAs are about 4.30 ± 0.07 dB below 
the corresponding QNL, which are shown in Fig. 3, a and b stand for quadrature-amplitude and quadrature-phase 
squeezing, respectively. 

Four NOPAs are locked individually by using Pound-Drever-Hall (PDH) method with a phase modulation of 56 
MHz on 1080 nm laser beam @. Each NOPA is operated at deamplification condition, which corresponds to lock 
the relative phase between the pump laser and the injected signal to (2n+l)7r (n is the integer). In the experiment, 
the relative phase locking is completed with a lock-in amplifier, where a signal at 15 kHz is modulated on the pump 
light by the piezo-electric transducer (PZT) mounted on a reflection mirror which is placed in the optical path of the 
pump laser and then the error signal is fed back to the other PZT which is mounted on a mirror placed in the optical 
path of the injected beam. 

In the beam-splitter network used to prepare eight-partite cluster states, Tl is phase-locked to zero, T2-T7 is 
phase-locked to 7r/2. The zero phase difference (Tl) between two interfered beams on a beam-splitter is locked by a 
lock- in amplifier with a modulation of 12 kHz. The 7r/2 phase difference (T2-T7) is locked by DC locking technique, 
where the photocurrent signal of the interference fringe is fed back to the PZT mounted on a mirror which is placed 
before the beam-splitter. 

In the homodyne detection system, zero phase difference for the measurement of quadrature-amplitude is locked by 
PDH technique with a phase modulation of 10.9 MHz on local oscillator beam. The 7r/2 phase for the measurement 
of quadrature-phase is locked by using DC locking technique too. 
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